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Abstract 

Previously, we have presented the first study to investigate the response 
of the Jovian thermosphere to transient variations in solar wind dynamic 
pressure, using a coupled, azimuthally symmetric global circulation model 



coupled with a simple magnetosphere model. This work (Yates et al. , 2012b) 
described the response of thermospheric flows, momentum sources, and the 
magnetosphere-ionosphere coupling currents to transient compressions and 
expansions in the magnetosphere. The present study describes the response 
of thermospheric heating, cooling and the auroral emissions to the aforemen- 
tioned transient events. We find that transient compressions and expansions, 
on time scales <3 hours, cause at least a factor of two increase in Joule heat- 
ing per unit volume. Ion drag significantly changes the kinetic energy of the 
thermospheric neutrals depending on whether the magnetosphere is com- 
pressed or expanded. These processes lead to local temperature variations 
>25 K and a ~2000 TW increase in the total power dissipated in the thermo- 
sphere. In terms of auroral processes, transient compressions increase main 
oval UV emission by a factor of ~4.5 whilst transient expansions increase 
this main emission by a more modest 37%. Both types of transient event 
cause shifts in the position of the main oval, of up to 1 ° latitude. 
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1. Introduction 



The Jovian upper atmospheric temperature is up to 700 K higher than 



that predicted by solar heating alone (Strobel and Smith, 1973 Yelle and 



Miller, 2004). This 'energy crisis' at Jupiter and the other giant planets has 



puzzled scientists for over 40 years. Different theories have been put forward 



to explain Jovian upper atmospheric heating: gravity waves (Young et al 



1997), auroral particle precipitation (Waite et al. , 1983; Grodent and Gerard 



2001), Joule heating (Waite et al. , 1983 Eviatar and Barbosa 1984) and ion 
drag QMiller et all 120001 ISmith et all 120051 iMillward et al.l 120051). None of 



the aforementioned studies have been able to fully account for the observa- 
tions. 



Magnetosphere-ionosphere (M-I) coupling models by Achilleos et al. ( 1998 ); 



Smith and Aylward (2009); Tao et al. (2009); Yates et al. (2012a) have all 



discussed steady-state heating and cooling terms in the Jovian thermosphere. 



Cowley et al. (2007) discussed 'transient' heating and dynamics in terms of 



power dissipated in the thermosphere via Joule heating and ion drag, as well 



as power used to accelerate magnetospheric plasma. Cowley et al. (2007) con 



sidered displacements of the Jovian magnetopause by ~40 Rj (Jovian radii; 
1 Rj=71492km). They found that, for compressions, there was a net trans- 
fer of power from magnetosphere to planet of ~325 TW, due to the expected 
super-corotation of magnetospheric plasma. For expansions, |Cowley et al. 



(2007) found that the power dissipated in the thermosphere (and used to 
accelerate magnetospheric plasma) increased by a factor of ~2.5 resulting 
from a large increase in azimuthal flow shear between the expanded magne- 
tosphere and the thermosphere. 



Melin et al. (2006) analysed infrared data from an auroral heating event 



observed by |Stallard et~aT] ( |2001[ |2002| > (from September 8-11, 1998) and 
found that particle precipitation could not account for the observed in- 
crease in ionospheric temperature (940— 1065 K). The combined estimate of 
Joule heating and ion drag increased from 67mWm~ 2 (on September 8) to 
277 mWm -2 (on September 11) resulting from a doubling of the ionospheric 
electric field (inferred from spectroscopic observations); this increase in heat- 
ing was able to account for the observed rise in temperature. Cooling terms 
(by Hydrocarbons and H% emission) also increased during the event but only 
by ~20 % of the total inferred increase in heating. Thus a net increase in 
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ionospheric temperature resulted. More detailed analysis showed that these 
cooling mechanisms would be unlikely to return the thermosphere to its ini- 



tial temperature before the onset of subsequent heating events. Melin et al 



(2006) thus concluded that the temperature increases could plausibly lead to 



an increase in equatorward winds, which transport thermal energy to lower 



latitudes (Waite et al, 1983) 



Cowley et al. (2007) modelled the auroral response of the Jovian system 



to transient compression and expansion events. 'Transient' in this context 
refers to events on timescales of <3 hours, during which the angular velocity 
of magnetospheric plasma is mainly determined by conservation of angular 
momentum. Cowley et al. (2007) found that the precipitating electron energy 



flux (~10 % of which is used to produce ultraviolet (UV) aurora) increases by 
two orders of magnitude at the open-closed field line boundary and decreases 
by ~50 % at the 'main auroral oval' during a transient magnetospheric com- 
pression. Essentially the opposite occurs for an expansion event, where main 
oval emission increases by a factor of 30 whilst emission at the open-closed 
boundary is reduced by ~98 %. Observations of transient aurorae have been 



carried out by Clarke et al. (2009) and Nichols et al. (2009). They showed 
that near the arrival of a solar wind shock at Jupiter's magnetopause (av- 
eraging an order of magnitude increase in solar wind dynamic pressure) the 



total UV auroral power doubled. Nichols et al. (2009) also showed that the 



location of the main oval shifted polewards by ~1 ° following the arrival of 
this shock, and that the increase in UV power was distributed evenly between 
the main oval and polar (high-latitude) emission. For solar wind rarefactions, 
these studies observed little or no change in total UV auroral power. 



Here, we use the Yates et al. (2012b) (henceforth Paper 1) model to 
simulate how the Jovian aurora and thermospheric heating responds to tran- 
sient variations in magnetospheric size (via changes in solar wind dynamic 
pressure). We employ this same methodology, as adapted in Paper 1; time- 
dependent plasma angular momentum profiles are used to represent changes 
in magnetospheric size. The response of thermospheric heating and aurora 
is subsequently investigated. The present work, along with Paper 1, are the 
first studies to model the Jovian thermospheric response to varying solar 
wind dynamic pressure using a global circulation model (GCM). 



In section [2] we summarise the scientific background appropriate for this 
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study, as well as model details. In sections [3] and [4] we present and discuss our 
findings for the transient compression and expansion scenarios respectively. 
We describe our conclusions in section [5] 



2. Scientific Background 

2.1. Time- dependence of the Jovian system 



Cowley and Bunce| ( [2003a] |b| and |Cowley et al.| Q2007[ ) were amongst the 



first studies to investigate the time-dependence of M-I coupling in the Jovian 
system. In this section we summarise the three time scales of relevance for 
this study which allow us to make simplifying assumptions. For a more de- 



tailed discussion the reader is referred Cowley et al. (2007) and Paper 1 



i) M-I coupling time scale: For the Jovian system the transfer of angular 
momentum from thermosphere to magnetosphere is estimated to occur on 
time scales of ~5— 20 hours by Cowley and Bunce (2003a). 



ii) Compression (and expansion) of the magnetosphere: Cowley and Bunce 



(2003a) and Cowley et al. (2007) estimated that magnetospheric size could 
change by as much as ~40 Rj in ~2— 3 hours. Since this figure is small com- 
pared to the M-I coupling time scale plasma angular momentum could be 
assumed constant during the modelled event. 

Hi) Thermospheric response time: There are currently two thermospheric 
response models proposed for transient changes in magnetospheric size: (i) 
Millward et all (120051) found, using a GCM flAchilleos et all 119981), that the 



thermosphere responds promptly (few tens of minutes), and (ii) thermosphere 



responds over a few days and is thus, considered to be unresponsive (Gong 



2005). We, however, make no assumptions about a thermospheric response 



time. We allow the GCM to respond independently to imposed variations in 
angular velocity of magnetospheric plasma. 

2.2. Magnetosphere model 

To investigate transient changes in magnetospheric size, we make use of 
two established magnetospheric models. These are: 

i) The quiescent magnetosphere model: This component is axisymmetric 



and consists of an equatorial magnetic field profile (Nichols and Cowley, 2004) 



and a magnetospheric plasma angular velocity profile, which is a steady-state 



4 



solution of the Hill-Pontius equation (Hill 1979 Pontius 1997). This qui- 



escent model has been presented in Nichols and Cowley (2004); Smith and 



Aylward (2009); Yates et al. (2012a) and Paper 1. As such, we will not dis- 



cuss it herein and we refer the reader to the cited studies for further details. 



ii) The transient magnetosphere model: This component is produced by 
perturbing the 'quiescent' equatorial magnetic field so as to predict its be- 



haviour during a transient compression or expansion. As described in Cowley 



et al. (2007) we use the the principle of magnetic flux conservation. We take 



the initial, quiescent field model and add a uniform southward (northward 
for expansion) perturbation field AB Z at each time step throughout the tran- 
sient event. AB Z is calculated, as a function of magnetodisc radius Rmm{P) 
(at each time step t) using the flux conservation condition as employed by 



Cowley et al. (2007) 



-7rR MM {t) 2 AB z = 2tt (F e {R MM {t)) - F e (R MM (t=0))). 



(1) 



Solving for AB Z gives 



AB * = JT%' (2) 

where F e is the initial equatorial flux function at time zero (i=0), AF = 
F e (RMM(t)) — F e (RMAi(t=0)) and AB z <0 for a southward field perturba- 
tion. t=0 denotes the initial, steady state of the system. 

Transient plasma angular velocity profiles ^m(-Rmm(^)) are obtained as- 
suming that plasma angular momentum is conserved. This is further dis- 
cussed in section 12.61 



2.3. Magnetosphere-ionosphere coupling 

The coupling of the magnetosphere and ionosphere models, has been dis- 



cussed in detail in previous studies by Smith and Aylward (2009); Yates et al 



(2012a) and Paper 1. Thus, here, we only briefly present equations relevant 
to the current study. For further details on the modelling of M-I coupling 
the reader is referred to the cited studies. 
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Electric and magnetic fields, along with ion-neutral collisions, create iono- 
spheric Pedersen currents flowing perpendicular to the direction of the plan- 
etary magnetic field at high latitudes (corresponding to the main auroral 
emission). These Pedersen currents are connected to the magnetosphere via 
field-aligned currents. Radial currents in the magnetosphere then close the 
'circuit' connecting the magnetosphere to the ionosphere. j\\i{0i) denotes 
field-aligned current (FAC) density just above the ionosphere at co-latitude 
9i. This quantity, when positive (upward from planet) corresponds to au- 
roral electron precipitation. Since FAC closes the divergence in horizontal 
Pedersen current, we may write: 



1 



dh 



J\\i 



2jcRf sin 9i dOi 



(3) 



The sign of j\\i{0i) indicates FAC direction (positive upward from planet), 
Ri is the ionospheric radius and Ip is the height-integrated Pedersen current 



(Cowley et al. 2007 Smith and Aylward 2009). 



The final aspect of M-I coupling we examine in this study is the energy 
transfer from planetary rotation to the thermosphere and magnetosphere. 
The energy transferred to the magnetosphere is used to accelerate magne- 
tospheric plasma towards corotation whilst the energy dissipated within the 
thermosphere is used for heating and increasing kinetic energy. The total 
power per unit area of the ionosphere transferred from planetary rotation P 
is the sum of atmospheric power Pa and magnetospheric power Pm dissipated 
per unit area ( Hill , 2001 ) . As shown by Smith et al. ( 2005 ) atmospheric power 
consists of two components: (i) Joule heating Pj and (ii) ion drag power P D , 
some of may be viscously dissipated as heat. These power relations are: 



p = 






Pm = 


VL m t, 




Pa = 






Pj = 






Pd = 




n T )r, 



(4) 
(5) 
(6) 
(7) 
(8) 



where 
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r = PiipBi 



(9) 



represents the torque exerted by the JxB force per unit area of the iono- 
sphere. Note that is a weighted average, computed over altitude, of ther- 
mospheric angular velocity (see Yates et al. (2012a) and Smith and Aylward 



(2009) for more detail). 



2.4- Auroral energies 

Once FAC densities have been calculated, 



we can use the methods of 



Knight (1973) and Lundin and Sandahl (1978), as presented in Cowley et al. 



(2007), to calculate the enhanced precipitating electron energy flux Ej 



E f 




+ 1 



(10) 



where Ef is the unaccelerated electron energy flux, j\\ i0 is the unaccelerated 
FAC density (or the maximum current that can be carried by the electrons in 
the absence of field-aligned potential drops) and is the upward (positive) 
FAC density calculated using Eq. ([3]). To enable a comparison with similar, 



earlier studies, we use the same electron population values described in Cow- 



ley et al. (2007), which are based on observations by Scudder et al. (1981) 



and Phillips et al. ( 1993a[b ). These parameters are presented in Table [T[ 



2.5. Thermospheric energy equations 

In this study we represent the thermosphere with a GCM. The GCM 
solves the non-linear Navier-Stokes equations of energy, momentum and con- 



tinuity, by explicit time integration ( Miiller- Wodarg et al. 



was originally created for Saturn (Miiller- Wodarg et al. 



Aylward, 2008) and later modified to represent Jupiter (Smith and Aylward 



2006). The GCM 



2006 Smith and 



2009). The model assumes azimuthal symmetry, and is thus two-dimensional 



(pressure/altitude and latitude). 



This study is focused on the response of atmospheric heating sources dur- 
ing transient events. Thus, we include only the atmospheric heating equa- 
tions in this section (Smith (2006): see also Achilleos et al. (1998): Tao et al. 



(2009)). For the Navier-Stokes equations on the thermospheric dynamics the 
reader is referred to Paper 1. The energy equation used in the GCM is: 



dip 
~dt 



q — u • V p ip — w 



dtp 
dp ' 



(11) 



with 



*P = T + U, (12) 
q = qp + qD + qj + q v + q K + q», (13) 

where ip is the sum of specific kinetic (T) and internal (U) energies of the 
system, u is the horizontal velocity and q is the sum of sources and sinks 
of heat. The second and third terms on the RHS represent advective energy 
terms. The individual components of q consist of the work done from: pres- 
sure gradients (qp), ion drag (?£>), Joule heating (qj), viscosity (q v ), thermal 
conduction (q K ) and solar radiation (q v ). w (w = dp/dt) is the convective 
time derivative of the atmospheric pressure p. The individual energy sources 
and sinks are expressed below: 



qp = 


_ w 
-u ■ gV p z H , 

Pm 




(14) 


qo = 


-u j x B, 

Pm 




(15) 


qj = 


— j • E*, 

Pm 




(16) 


q v = 


2 d ( du\ 




(17) 


q K = 




dT\ 
dp J ' 


(18) 



where g is the acceleration due to gravity, z is altitude (above 300 km), p m is 
the thermospheric density, j, B and E* are respectively the current, magnetic 
field strength and rest-frame electric field. Finally, rj is the viscosity, k is 
the thermal conductivity and T the temperature. The absorption of solar 
radiation is a small energy source in the current model (Smith, 2006). For 
details regarding the absorption of solar radiation the reader is referred to 



Moore et al. (2004) and Miiller-Wodarg et al. (2006). 
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2.6. Obtaining the transient plasma angular velocity 

As the calculation of plasma angular velocity in the transient regime has 
been discussed thoroughly in Paper 1, we shall only briefly summarise the 



main points in this section. The reader is referred to Paper 1 and/or Cowley 



et al. (2007) (on which our model is based) for a more detailed description 



of the transient magnetospheric angular velocity model. 

We may assume conservation of the plasma's angular momentum if we 
consider a solar wind shock or rarefaction which will cause the magnetosphere 



to compress or expand several tens of Rj on time scales of ~2— 3 hours (Cow- 



ley et al. (2007) and Paper 1). This gives a plasma angular velocity profile 



^Af(#j,£), at any time t (throughout the 'pulse' in solar wind pressure) of 



n*(M) = n M (M=o) ( Pe{6 ! , t ift Y, (w) 

V Pe{t>i,t) J 

where the notation t=0 and t denotes the initial (steady-state) and transient 
(at each time-step throughout the event) states respectively. p e (0i,t) is the 
time-dependent equatorial radial distance which magnetically maps to a co- 
latitude 9i. Note that the size of the magnetodisc Rmm{P) in this study is 
determined using a Gaussian function (see section 2.6 in Paper 1) and that 
the full duration of the solar wind pulse is three hours. In sections [3] and 
[4] we present results pertaining to the initial steady-state, half-way through 
the pulse (point of maximum compression or expansion) and at the end of 
the pulse (where -Rmm(^) returns to the initial value .Rmm(£=0)). 



The magnetosphere in our model is divided into four regions (Yates et al. 



(2012a) and Paper 1). Region I represents the open field lines mapping to the 
polar cap; region II represents the outer magnetosphere; region III represents 
the middle magnetosphere (magnetodisc); and region IV is the inner magne- 
tosphere, which we assume is fully corotating in steady state. Region III is 
shaded light grey as this is our main region of interest due to its dominant 
role in M-I coupling. 

Plasma angular velocities (thin lines) and corresponding thermospheric 
angular velocities (thick lines) as functions of ionospheric latitude are shown 
in Fig. [TJ This Figure and the system state labels are identical to Figure 2 
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from Paper 1; thus, we will only summarise the transient plasma angular 
velocity discussion. 



Angular velocity profiles for our transient compression event are shown 
in Fig. [l^,. The solid lines indicate our initial steady state, henceforth called 
'case CS' (pre-Compression Steady-state). Dashed lines represent results 
midway through the event (-Rmm(^) is at a minimum) and are referred to as 
case CH (Compression Half-pulse). The dot-dashed lines indicate the sys- 
tem at the end of the transient pulse, case CF (Compression Full-pulse). 
Considering case CH, fl M super-corotates over most of regions IV and III. 
This has a number of consequences such as the corresponding increase in 
Qt, ionospheric current reversal in these regions and the reversal of momen- 



tum transfer between the thermosphere and magnetosphere (Cowley et al. 



2007). The magnetospheric response from case CS to CH is rather large 



(factor of two increase in peak Q M ) compared to the ~3 % average increase 
in flrp] this is due to the neutral thermosphere having a larger mass than the 
magnetospheric plasma region to which it maps, and thus a larger inertia. 
For case CF, n u (CF)=n M (CS) but tt T {C F)^Sl T {C S); this demonstrates 
that there is a time lag in the response of the thermosphere compared to 
the magnetosphere during transient changes in solar wind dynamic pressure. 
After the pulse has subsided, the thermosphere remains perturbed relative 
to its steady state. 



Fig. [T]d shows angular velocity profiles corresponding to the transient 
expansion scenario. System state labels and line styles are as in Fig. [1^ but 
the 'C (Compression) is exchanged for an 'E' (Expansion) giving cases ES, 
EH and EF. In case EH, magnetospheric sub-corotation decreased drastically 
which consequently decreased Qt, albeit to a much lesser degree. As in the 
compression scenario, Qm(EF)=Qm(ES) and Qt(EF)^Qt(ES). We do, 
however, note that a small amount of current reversal is expected near the 
region III/II boundary due to VL M (EF)>Vt T (EF)] this acts to effectively 
'pull' the thermosphere back towards corotation after a period of substantial 
sub-corotation. 



3. Magnetospheric Compressions 

In this section we present auroral and thermospheric heating results for 
a transient magnetospheric compression event. We present results from 
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three phases throughout the transient event: i) Steady state (case CS with 
-Rmm=85Rj), ii) 'half-pulse' (case CH with _Rmm=45Rj) and, iii) 'full-pulse' 
(case CF with i?MM=85Rj). 

3.1. Auroral response 

We begin by discussing how a transient magnetospheric compression af- 
fects the aurorae. Here we only discuss the effects on FAC densities and the 
energy flux carried by precipitating electrons. For the detailed response of 
global M-I coupling currents the reader is referred to Paper 1. 



We begin our discussion by briefly summarising the response of FAC 
densities (see Paper 1 for detailed derivation). Fig. [2^ shows FAC density 
as a function of latitude for cases CS (solid line), CH (dashed line) and CF 
(dot-dashed line, nearly coincident with CS). Magnetospheric regions (region 
III shaded) are separated by black dotted lines and the latitudinal size of a 
Hubble Space Telescope (HST) ACS-SBC pixel (0.03x0.03 arc sec) is repre- 
sented by the dark grey rectangle (assuming that the magnetic axis of the 
Jovian dipole is perpendicular to the observer's line of sight). We notice 
that between cases CS and CF there are only minor differences in FAC pro- 
files. These differences result from the lag in response of the thermosphere to 
the transient pulse. For case CF, tt M {CF)=tt M {CS) but n T (CF)^n T (CS), 
meaning the thermosphere has not settled back to its steady state profile (due 
to its large inertia; see section 2.6 and Paper 1). Comparing case CS with 
CH, CH has: i) substantial current reversal in region III, ii) maximum FAC 
about twice that of case CS, iii) increased FAC at the region II/I boundary 
(although this result is sensitive to our assumed choice of uniform angular 
velocity in region II). The current reversal over ~4° of region III is a direct 
result of the super-corotation of Hm (see Fig. [T^,) whilst the increases in up- 
ward FAC result from increases in the Pedersen current (see Paper 1) caused 
by an increase in flow shear (Q T — Qm)- 



Corresponding precipitating electron energy fluxes are shown in Fig. [2}d. 
These fluxes are plotted as a function of latitude. The line style code and 
labels are the same as in Fig. [2^i and the grey solid line indicates the limit 
of present detectability with HST instrumentation 
( |20"07| )). These fluxes are obtained using Eq 
ward (positive) FAC densities presented in Fig 



'1 kR; Cowley et al. 



(10) 



which uses only the up- 
2k. We initially compare 
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electron energy fluxes for cases CS and CF. These profiles are identical pole- 
ward of ~74°; equatorward of this location, case CF shows little deviation 



from CS, despite the thermospheric lag discussed in section |2.6| In region 
III, we find that the peak energy flux for case CF is ~35 % larger than that 
in case CS and the location of these peaks coincides with the location of the 
main auroral oval (~74°). The slight increase in peak energy flux is due to 
a relative increase in flow shear as seen in Fig. [T^i. Case CF would therefore 
produce main oval emission of approximately ~200 kR brighter than that of 
CS (assuming that lmWm" 2 of precipitation creates ~10kR of UV output 



(Cowley et al. 2007)). The main oval would also be ~1.5° broader for CF 



than for case CS. 

The Ef profile for case CH is different from those of both cases CS and 
CF. There are three main changes in CH compared to CS: i) peak energy 
flux in region III is ~280mWm" 2 , almost a factor of five larger, ii) location 
of peak energy flux has shifted polewards by ~0.2° and, iii) presence of a 
second peak with an energy flux of 1.7mWm~ 2 at the region II /I bound- 
ary. The large increase in energy flux is caused by a substantial increase in 
flow shear between the thermosphere and magnetosphere, resulting from the 
super-corotation of the middle magnetosphere (see Fig. [TJ. The presence of 
a second peak at the region II/I boundary is also due to flow shear increase 
across the boundary, as the magnetosphere in region II corotates at a larger 
fraction of Qj compared to case CS. The result for this higher-latitude bound- 
ary should be regarded as preliminary, since it is sensitive to the values of Qm 
we assume in the outer magnetospheric region and polar cap. Flow velocities 



in these regions are poorly constrained by observations (Stallard et al. , 2003) 



Comparing the energy flux profile of case CH with the equivalent case in 



Cowley et al. (2007) (with a responsive thermosphere; they do not calculate 
energy fluxes for a non-responsive thermosphere), we see that in the closed 
field regions (III and II), peak energy fluxes are two orders of magnitude 
larger in case CH. This demonstrates the differences between using a GCM 
to represent the thermosphere or using a simple 'slippage' relation between 
thermospheric and magnetospheric angular velocities. At the open-closed 
field line boundary (II/I boundary) our peak flux is an order of magnitude 



smaller than that in Cowley et al. (2007); this difference arises from the 



different models used to represent the outer magnetosphere. The outer mag- 
netosphere (region II) and open field line region (region I) in this study is 
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modelled using Cowley et al. (2005). 



The increase in Ef for case CH would lead to corresponding increases in 
auroral emission. As such, we would expect 'main oval' emission for case CH 
to shift polewards by ~0.2 ° and be ~4.7 x larger than emission in case CS i.e. 
~2800 kR compared to ~600 kR. IClarke et all (120091) andlNichols et al.l ([20091 



observed that the brightness of UV auroral emission increased by a factor of 
approximately two, in response to transient (almost instantaneous) increases 
in solar wind dynamic pressure (~0.01— 0.3 nPa or equivalently ~109— 72 Rj). 



Nichols et al. (2009) also observed poleward shifts in main oval emission on 
the order of ~1 ° corresponding to the arrival of solar wind shocks. Total 
emitted UV power can be used to describe auroral activity, assuming that 
this quantity is ~10% of the integrated electron energy flux per hemisphere 
(Cowley et al. , 2007). Case CH has a total UV power of ~ 1.58 TW (compared 



to ~420 GW for case CS), which is a factor of two to three times larger 



than UV powers observed by both Clarke et al. (2009) and Nichols et al 



(2009). The profile of case CH also indicates the possibility of observable 



polar emission at region II /I (open-closed) boundary. This conclusion is, 
however, sensitive to our model assumptions. 

3.2. Thermospheric heating 

In this section, we present and discuss the response of thermospheric heat- 
ing terms to our transient magnetospheric compression event. We present 
azimuthal and meridional flows and temperature differences (from steady 
state) in Fig. [3] The results in Fig. [3] have been discussed at length in Paper 
1 and we therefore only describe the major flow patterns associated with 
zonal (Figs. |3^i-c) and meridional flow (Figs. [3]i-f). Thermospheric tempera- 



i and we use thermospheric heating terms and 
4) to aid in their discussion. For completion we 
also present integrated powers per hemisphere in Fig. [5] Cases CS-CF are 
represented in columns 1—3 of Figs. [3]j4] respectively. 



tures are presented in Figs, 
powers per unit area (Fig. 



Figs. [3^i-c show thermospheric azimuthal velocity, in the corotating ref- 
erence frame, as a function of altitude and latitude for cases CS-CF. Pos- 
itive velocities indicate super-corotation whilst negative velocities indicate 
sub-corotation. The white contour indicates the locus of rigid corotation, 
black contour indicates regions of strong super-corotation (> 25ms -1 ), the 
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white dashed contour shows regions that are strongly (> — 2500 m s -1 ) sub- 
corotating and the black arrows indicate the direction of meridional flow. 
Magnetospheric regions are labelled and separated by black dotted lines. 
The zonal flow patterns include: i) a low-altitude (~400 — 500 km) super- 
corotating jet, centred at ~72°; and ii) a large sub-corotating jet, spanning 
region III - I (blue region in Figs. [3^-c). These general features are persistent 
throughout the transient compression event. Differences do arise in the speed 
of flows: where we see an increase in the degree of corotation (and super- 
corotation) in case CH (indicated by the area covered by the white-dashed 
contour having decreased and that covered by the black contour having in- 
creased). This is followed by a subsequent decrease super-corotation in case 
CF. We do however note that in regions II, there is a slight increase in coro- 
tation in case CF. Figs. [3]i-f show thermospheric meridional velocity as a 
function of altitude and latitude for cases CS-CF. The arrows indicate the 
direction of the steady-state meridional flow pattern (Fig.[3]i), where we have 
low altitude (<600km) poleward flows and high altitude (>600km) equator- 
ward flows, as predicted by Smith et al. (2007); Smith and Aylward (2009) 
and Yates et al. (2012a). We also note the presence of localised accelerated re- 



gions (red/brown hues) which result from the advection of momentum (Yates 



et al. , 2012a). The response of the meridional flows can be seen in Figs. 3tl-f 



where flow patterns evolve with time to become quite complex. The number 
of local regions with strongly accelerated flows increases significantly with 
time. The reader is referred to Paper 1 for more detailed discussions on 
these flows and their causes. 



Fig. [3J5 shows thermospheric temperature as a function of altitude and 
latitude for case CS. Figs. [3ji-i show the difference in thermospheric tempera- 
ture between cases CH and CS, and cases CF and CS. We will use Figs, [pi-f, 
showing contour plots for various thermospheric heating (Figs. |4^l-c) and 
cooling (Figs. |4]i-f) terms (see plot legends for details) to interpret the tem- 
perature response. 



In Fig. [3J5 we see a clear temperature difference between upper (>75°) 
and lower (<75°) latitudes; lower latitudes are cooled whilst upper latitudes 
are significantly heated ( |Smith et al. 2007; Smith and Aylward 2009). We 
see a 'hot spot' (in region I) with a peak temperature of ~705K. This arises 
from the poleward transport of Joule heating (from regions III and II) by 
the accelerated meridional flows shown in Fig. [3]i (Smith and Aylward |2009 
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Yates et al. 



2012a). 



Fig. [3)2 shows the temperature difference between cases CH and CS. There 
are three prominent features in Fig. |3ji: 

i) Temperature increase up to ~26K across the region III/II boundary 
(^>400km) resulting from a large (x2) increase in Joule heating and the 
addition of other heat sources, such as adiabatic heating (see Fig. |4Jd) . The 
large increase in Joule heating is caused by the increase in FAC and the 
rest-frame electric field (see Fig. [2^, and Paper 1 respectively). 

ii) Temperature decrease down to ~— 22 K, at low altitudes of region 
II. Fig. |4}d shows that at low altitudes (<500km) of region II there is, on 
average, a 20 % decrease in energy deposited by Joule heating and ion drag. 
This, coupled with the presence of energy lost by ion drag (Fig. |4^) in this 
region causes the significant decrease in temperature shown above. All the 
factors discussed above result from the reversal and decrease (in magnitude) 
of the flow shear between the magnetosphere and thermosphere in case CH. 

Hi) A maximum of ~17K increase at low altitudes in region I. The merid- 
ional velocity of case CH increases slightly (~2 %) in in this region and, as 
such, can transport heat from Joule heating and ion drag polewards more 
efficiently than in case CS. 



Fig. [3} shows the temperature difference between cases CF and CS. Im- 
mediately, we can see that there are significant changes in the distribution of 
temperature in the upper thermosphere of case CF. There are four 'finger- 
like' regions with local temperature increases >50 K (maximum of 175 K) and 
three regions with temperature decreases <40 K. These alternating temper- 
ature deviations increase with altitude and are collocated with accelerated 
meridional flow regions. Considering Figs. [4^ and f, we see that the heat- 
ing and cooling terms are now quite complex, with advective and adiabatic 
terms dominating (>10 x Joule heating and ion drag terms). The CF ther- 
mosphere appears to be transporting heat (in 'waves'), both equatorward and 
poleward from the region III/II boundary (see Fig.[3]f). Achilleos et al. (1998) 



also shows a similar phenomenon (see top left of Fig. 9 in Achilleos et al. 



( p98) )), whereby ' wave-like' patterns of high temperature are being trans- 
ported away from the auroral region. The energy deposited in the auroral 
regions heats the local thermosphere which increases local pressure gradients. 
Advection then attempts to redistribute this heat which momentarily cools 
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the local area until enough heat is deposited again and the process restarts. 



Figs. |4^;-i show powers per unit area as functions of ionospheric latitude, 
calculated using Eqs. ([5] -[8]), for cases CS, CH and CF, respectively. In- 
tegrating these powers over the hemisphere gives the global powers shown 
in Fig. |5j Powers are indicated by legends on figures. In case CS, magneto- 
spheric power (power used to accelerate magnetospheric plasma) is dominant 
up to ~73°, where atmospheric power (sum of Joule heating and ion drag) 
quickly dominates for all poleward latitudes (see Fig. [IJ2;). This indicates 
that relatively expanded M-I systems (in steady-state) generally dissipate 
more heat in the atmosphere than in acceleration of outward-moving plasma 



(Yates et al. 2012a). For case CH, Fig. [4ji, we see the effects of plasma 
super-corotation in region III, where magnetospheric power reverses (now 
negative), heat dissipated as Joule heating doubles (as shown above), pos- 
itive ion drag decreases by ~70 % and negative ion drag increases by two 
orders of magnitude. Powers decrease in region II due to the decrease in 
flow shear (see Fig. [T^i). In case CF, powers per unit area closely resemble 
those for case CS. There are small increases in peak magnetospheric power 
(~10 %) and Joule heating (~25 %) leading to an overall maximum increase 
in available power of ~10 %. These findings are represented on a global scale 
in Fig. [5j 

Overall, midway through a transient compression, magnetospheric plasma 
super-corotates, which leads to energy being transferred from the magneto- 
sphere to the thermosphere. Our results indicate that ~2000TW of total 
power (magnetospheric and atmospheric) is gained by the coupled system 
as a result of plasma super-corotation. Note that this is considerably larger 



than the ~325 TW (closed and open field regions) calculated in Cowley et al. 



(2007) for a responsive thermosphere scenario. This energy transfer from the 



magnetosphere would act to, essentially 'spin up' the planet (Cowley et al. 



2007) and increase the thermospheric temperature. In case CF, plasma is not 
super-corotating; thus the picture is fairly similar to case CS. The main dif- 
ference is that there is a ~20 % increase in total power dissipated in the atmo- 
sphere and in acceleration of the magnetosphere. This arises from increases 
in flow shear due to the 'lagging' thermosphere (see Fig. [T]) and inevitably 
leads to the local temperature increases seen in Fig. [3} and discussed above. 
The finite time required for thermospheric response results in the described 
'residual' perturbations to the initial system (CS) even after the pulse has 
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subsided (CF). 



4. Magnetospheric Expansions 

In this section we present auroral and thermospheric heating results for 
a transient expansion event. We present results from three periods through- 
out the transient event: i) Steady state (case ES with R MM =45 Rj), ii) 
'half-pulse' (case EH with Rmm= 85 Rj) and, iii) 'full-pulse' (case EF with 
-Rmm=45 Rj). 

4-1. Auroral response 

We begin by investigating how a transient magnetospheric expansion af- 
fects the aurorae. Thus, we only discuss the effects on FAC densities and the 
precipitating energy flux. For the detailed response of M-I coupling currents, 
the reader is referred to Paper 1. 

As for the case of our compression scenario, we summarise here the re- 
sponse of FAC densities (see Paper 1 for detailed discussion). Fig. [6^1 shows 
the variation of FAC density with ionospheric latitude for cases ES (solid 
line), EH (dashed line) and EF (dot-dashed line). Magnetospheric regions 
(region III shaded) are separated by black dotted lines and the latitudinal 
size of a HST ACS-SBC pixel is indicated by the dark grey rectangle. Case 
EF exhibits only minor differences with case ES, which result from the lag 
in the thermospheric response to transient magnetospheric changes. For case 
EH, the magnetospheric expansion causes the magnetic field strength (in the 
magnetodisc) and plasma angular velocity to decrease. Flow shear (f2y— Om) 
thus increases and results in a broad region of upward FAC from ~69 to 
~74° latitude. We also note a ~15% increase in peak FAC density in this 
broad upward region. At the region III /II boundary there is a small region of 
strong downward FAC produced by large gradients in Qm and Ep across this 
boundary. We also note that poleward of region III, only downward FACs 
exist due to the small imposed change in Qm in regions II and I. 



We interpret our FAC density profiles by considering the corresponding 
precipitating electron energy fluxes, shown in Fig. [6Jd. Fluxes are plotted 
as functions of latitude. The line style code and labels are the same as in 
Fig. [6^ and the grey solid line indicates the limit of present HST detectability 

(2007)). We begin by comparing cases ES and EF, 



1 kR; Cowley et al. 
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which are almost identical except between ~71 ° and ~73 ° latitude. These 
differences are caused by the finite thermospheric response time. Both pro- 
files peak at ~74° latitude, equivalent to the location of the 'main auroral 
oval', and at ~80°, the boundary between open (region I) and closed field 
lines (region II). Therefore, we would expect a fairly bright auroral oval of 
~88kR for case ES and ~79kR for case EF. The electron energy flux for 
case EF is ~10 % smaller than case ES due to Qt(ES)>Qt(EF) leading to 
a smaller flow shear. Our model also predicts the possibility of observable 
polar emission (region II /I boundary) of ~15kR for both cases ES and EF. 
However, this region is strongly dependent on the model used and poorly 
constrained by observations. 



We now discuss the response of the energy flux Ef midway through our 
transient expansion event. The case EH profile, poleward of ~74° latitude, 
lies entirely below the limit of detectability due to the downward (negative) 
FAC density in this region. In region III, there are now two peaks, separated 
by ~1°. The first peak, located at ~73° is ~37% larger than the second 
one, at ~74°. These peaks result from the large degree of magnetospheric 
sub-corotation in region III (evid ent in Fig. [T]o) . Comparing case EH with 
the equivalent expansion case in Cowley et al. (2007), case EH, in region 
III, has a maximum value of Ef that is twice that in Cowley et al. (2007). 
This study represents the thermosphere with a GCM which responds inde- 
pendently to VLm profiles. Our results indicate that this response is not as 
strong as that in Cowley et al. (2007), who use a simple 'slippage' relation 
to model the thermospheric angular velocity. At the open-closed field line 
(region II/I) boundary, Cowley et al. (2007) obtain large energy fluxes due 
to their large change in VLm across these regions; in our study, we obtain 
negligible changes in Ef due to our smaller change in imposed Qm across 
this boundary. 



Considering the auroral response in case EH, we expect to see no polar 
emission (a darkening from case ES) at the region II/I boundary. In region 
III, we expect a ~1 ° equatorward shift along with a small increase (~40kR) 
in peak main oval brightness compared to case ES. We also note the possible 
observation of a somewhat bifurcated main oval; with emission peaking at 
~73° and ~74° latitude. The main oval would also be considerably broader 
(~2— 3°) as a result of the large increase in the spatial region of magneto- 
spheric sub-corotation. Clarke et al. (2009) observed little change in auroral 



18 



brightness near the arrival of a solar wind rarefaction region, however Nichols 



et al. (2009) have seen changes in main oval location. The total UV power in 
case EH is ~270 GW (compared to ~78 GW in case ES). While this power 
is considerably smaller than that in case CH, it is comparable to UV powers 



calculated in Clarke et al. (2009) and Nichols et al. (2009), following solar 



wind rarefactions (~200 - 400 GW).. 

4-2. Thermospheric heating 

In this section, we present and discuss the response of thermospheric 
heating terms to our transient expansion event. We present azimuthal and 
meridional flows and temperature differences in Fig. [7j The results in Fig. [7] 
have been discussed at length in Paper 1, and we therefore only describe 
the major flow patterns associated with zonal (Figs. (7^,-c) and meridional 
directions (Figs. [7}d-f). Thermospheric heating terms and powers per unit 
area are presented in Fig. [8] and we use Figs. [7^;-i to aid in their discussion. 
For completeness, we also present integrated powers per hemisphere in Fig. [9] 
Cases ES, EH and EF are represented in respective columns 1—3 of Figs.[7jj8} 



Figs. [7^i-c show the variation of thermospheric azimuthal velocity, in the 
corotating reference frame, with altitude and latitude for cases ES, EH and 
EF. Positive velocities represent super-corotation whilst negative velocities 
indicate sub-corotation. The white contour indicates the locus of rigid coro- 
tation and the black arrows indicate the direction of meridional flow. Mag- 
netospheric regions are labelled and separated by black dotted lines. The 



zonal flow patterns are as in section 3.2 i.e a low altitude super-corotating 
jet in region III and a large sub-corotating jet spanning regions III-I. We 
note that the area experiencing strong super-corotation (black contour) is 
at a minimum in case EH and that the degree of sub-corotation decreases 
as the system responds to the event. Figs. [7]d-f show thermospheric merid- 
ional velocity as a function of altitude and latitude for cases ES-EF. Labels, 
arrows and contours are the same as in Figs. [7^i-c. The steady-state flow pat- 
tern, consisting of high-altitude equatorward flows and low-altitude poleward 
flows, is as in case CS (discussed in Smith et al. ( |2007 ); Smith and Aylward 
(2009); Yates et al. (2012a) and Paper 1). As the system progresses through 



the expansion event, the high altitude meridional flow reverses for parts of 
regions III and II and all of I. Large, localised acceleration regions also de- 
velop. For further details on thermospheric flows, the reader is referred to 
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Paper 1. 



Fig. [7^; shows temperature as a function of altitude and latitude for case 
ES. Figs. [7ji-i show the difference in thermospheric temperature between 
cases EH and ES, and cases EF and ES, as functions of altitude and latitude. 
Magnetospheric regions are labelled and separated by black dotted lines and 
temperatures are indicated by the colour bar. We interpret the response of 
thermospheric temperature with the aid of Figs. [8^i-f, showing contour plots 
for various thermospheric heating (Figs. [8^,-c) and cooling (Figs. [8jl-f) terms 
(see plot legends for details). 



Fig. [7J5 shows similar results to those described in section 3J2 The main 



difference is related to the polar 'hot spot' which is considerably cooler (peak 
temperature of ~590K) than that for case CS (peak temperature ~705K). 
As previously discussed, the 'hotspot' results from the meridional trans- 
port, via poleward accelerated flows, of Joule heating from lower latitudes 
(~73- 84°; see Figs, and d) flSmith and Aylward[ |2009[ |Yates et al. 



2012a). 



Fig. [7ji exhibits the temperature difference between cases EH and ES. The 
figure shows a ~50K local temperature increase at low altitudes (<700km) 
in regions III and II. This is caused by a >4x increase in Joule heating and 
ion drag power in this region (see Figs. [7]d and e). Also evident are two more 
minor temperature variations: i) ~10K decrease at high altitude, centred 
on the region III /II boundary and, ii) ~10K increase in the polar 'hotspot' 
region. We interpret that this local 'cold spot' arises mainly from the appear- 
ance of meridional accelerated flows, which efficiently transport heat away 
(equatorward and poleward) from the region III/II boundary, and the factor 
of three increase in adiabatic cooling. We attribute the small increase in tem- 
perature at the polar 'hotspot' to the 35 % increase in low altitude meridional 
velocity (itself due to poleward acceleration - more detail in Paper 1). 



Fig. [7} shows the temperature difference between cases EF and ES. The 
temperature profile has changed significantly from that in Fig. (7}i. There are 
two regions where temperatures increase by up to ~50 K: i) extending from 
~73— 85 latitude and low altitudes in regions III and II, and all altitudes in 
region I (these map to the large poleward-accelerated region in Fig. [7]f); and 
ii) high-altitude (>600km) region, centred at ~66° latitude. These regions 
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are primarily heated by horizontal advection (high-altitude only) and adia- 
batic terms (all altitudes) as shown in Fig. [8J3; these terms have increased 
from case ES by, at most, 800 % and 500% respectively. The final feature of 
note in Fig. [7| is the region cooled by up to ~— 22 K, lying between the two 
heated regions at altitudes >550 km. This cooling is caused by a combination 
of local increases in horizontal advection and adiabatic cooling, by factors of 
three and greater. Similar to case CF, case EF's meridional flows seem to be 
transporting heat equatorward and poleward, although the majority of these 
flows act to transport thermal energy poleward. 



Figs. |8^-i show powers per unit £1X621 cLS functions of ionospheric latitude 
for cases ES-EF respectively. Integrating these powers over the hemisphere 
gives the global powers in Fig. |9j Power sources are indicated by the legends 
on the figures. Fig. [8J5 shows the energy balance in the thermosphere for case 



ES. As discussed in Yates et al. (2012a), most of the energy in region III is 



expended in accelerating magnetospheric plasma; in region II we have a situ- 
ation where magnetospheric power and atmospheric power (the sum of Joule 
heating and ion drag) are equal, due to Qm=0.5Qj. Atmospheric power is 
dominant in region I. For case EH (Fig. [8|i), the magnetospheric power is 
dominant up to ~73 latitude; atmospheric power, primarily Joule heating, 
then dominates at all poleward latitudes. This results from the large increase 
in flow shear (Q? — ^m)- The magnetosphere of case EF super-corotates, 
compared to the thermosphere, at latitudes <73°. This causes a reversal 
in energy transfer, which now flows from magnetosphere to atmosphere (see 
Fig. [8}). Polewards of 73°, the energy balance is similar to that of case ES. 

We conclude that a transient magnetospheric expansion creates a signifi- 
cant increase in both power dissipated in the atmosphere due to Joule heating 
(~6x that of ES) and ion drag (~3x that of ES). Moreover, the power used 
to accelerate the magnetosphere towards corotation is ~7x that of ES, and 
is shown in Fig. [9] These increases lead to a total power per hemisphere 
of ~2600 TW which is three times larger than the responsive thermosphere 



case in |Cowley et al. (2007). These changes in heating and cooling create the 
local temperature increases discussed above, but only by <7K (compared 
to case ES). For case EF, where we now have the magnetosphere rotating 
faster than the thermosphere, there is a ~75 % decrease in the magnitude 
of 'magnetospheric' power. The magnetosphere is thus transferring power to 
the thermosphere in this case, albeit a relatively small amount. This effec- 
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tively 'pulls' the thermosphere along, increasing its angular velocity in order 
to return to the steady-state situation where Qt>^m- We note that energy 
dissipated via Joule heating also decreases slightly due to the small decrease 
in flow shear. Overall, then, the total power per hemisphere in case EF is 
only 30 % that of case ES. 



5. Conclusion 

This study models the effects of transient changes in solar wind dynamic 
pressure on Jovian auroral parameters and the thermospheric energy bud- 
get. A companion study, considering the response of M-I coupling currents 
and thermospheric dynamics is presented in Paper 1. The transient changes 
simulated are two types: i) a transient compression event, and ii) a transient 
expansion event. Both events are assumed to last for three hours, and we 
have presented and compared pre-event, midway (maximum compression or 
expansion) and post-event results in sections [3] and El 



Simulations of the auroral response to a transient compression and the 
subsequent return to the initial magnetospheric configuration showed sig- 
nificant variation in auroral size, location and brightness. Midway through 
the compression event we find that auroral brightness, estimated from the 
precipitating electron energy flux, increases from ~600 kR (in steady state) 
to ~2800 kR. The location of the main auroral oval also shifts poleward by 
~0.2° latitude and its width shrinks slightly compared to the steady-state. 
Observations by Clarke et al. (2009) and Nichols et al. (2009) find factor-of- 
two increases in auroral brightness as well as poleward shifts of up to ~1 ° 
(Nichols et al. , 2009) due to the arrival of solar wind shocks. The total UV 
power can also be used to examine the auroral response, and we find that 
UV powers are approximately a factor of two or three larger than those ob- 
served by Clarke et al. (2009) and |Nichols et aL (2009). This is interesting 
as their observed compressions (magnetopause change of ~40Rj) are larger 
than our modelled one (magnetopause change of ~33Rj). Immediately af- 
ter the transient compression event, we find that the aurora differs slightly 
from the initial steady-state, and then only in regions mapping to the middle 
magnetosphere. There is an increase in auroral emission of ~200 kR (~35 %) 
and a ~1.5° broadening of the main oval. 
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The super-corotation of magnetospheric plasma midway through the com- 
pression event caused various changes in thermospheric heating and cooling 
terms. Joule heating increases by a factor of two and negative ion drag (de- 
creases kinetic energy of the system) increases by a factor of seven, whilst pos- 
itive ion drag (increase kinetic energy of the system) is halved. Overall, the 
magnetospheric super-corotation results in ~2000 TW resultant power being 
dissipated in the thermosphere. This leads to local temperature changes of 
~25 K. After the transient event subsides, adiabatic and horizontal advective 
terms dominate over Joule heating and ion drag by an order of magnitude. 
This results in high-altitude thermospheric temperature increases of >40 K. 



Our transient expansion simulations show less auroral variation than the 
compression simulations. At maximum expansion, peak auroral brightness 
increases from ~88kR to ~126kR and shifts equatorwards by ~1 ° latitude. 
We do note ,however, the presence of a slightly smaller auroral peak, ap- 
proximately 1 ° poleward, creating a slight bifurcation of the main oval. The 
latitudinal width of the main oval also increased by ~2— 3°. We would thus 
expect to observe a broad, slightly 'bifurcated' main oval and no significant 
polar emission midway through the expansion event. Observations of the 
auroral response near solar wind rarefaction regions by Clarke et al. (2009) 
have shown that there is little, if any, corresponding change in auroral bright- 
ness. As the expansion event subsides, we see only a small difference between 
the initial auroral profile and the final one. This difference is a small 10 % 
decrease in peak auroral emission, caused by the 'lagging' thermospheric re- 
sponse. 



The transient magnetospheric expansion, similar to the compression, in- 
creased the total energy available to the thermosphere-magnetosphere sys- 
tem. The expansion lead to a factor-of-four increase in ion drag power and 
Joule heating rate (in the region mapped to the middle magnetosphere), 
which resulted in local temperature increases <50 K. Globally, the large 
degree of magnetospheric sub-corotation increased the power expended in 
accelerating the magnetospheric plasma by a factor of seven. This aspect 
coupled with the changes in Joule heating, produced a ~2170TW increase 
in the total power dissipated and used by the thermosphere. At the end of 
the expansion event, horizontal advection and adiabatic heating terms domi- 
nate, and are on average 650 % larger than in steady-state. This creates large 
areas of temperature variation <50 K. Globally, due to the magnetosphere 
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rotating slightly faster than the thermosphere, the total power dissipated in 
the coupled system decreases significantly to ~125TW (~30% of the total 
in steady-state). 



In summary, by examining our model's meridional flows and temperature 
distributions, we find that, by the end of the imposed transient changes, heat 
(thermal energy) has started to be transported both poleward and equator- 



ward, particularly at high altitudes, in accordance with studies by Waite 



et al. (1983) and Melin et al. (2006). At present, though promising, these 



results are unable to explain the high thermospheric temperatures observed 
at low latitudes in Jupiter (~900K; Seiff et al. (1998)). As part of future 
studies, it would be interesting to vary the duration of these transient events 
and to experiment with a more realistic solar wind profile, in order to model 
the dynamic response of the Jovian thermosphere over more extended periods 
of time. 
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Table 1: Magnetospheric electron source parameters. This table is adapted from Table 1, 
Cowley et all (120071). N e represents the electron density, W t h the electron thermal energy, 



the unaccelerated current density and Ef the unaccelerated energy flux. 



Parameter 


Open field lines 


Outer magnetosphere 


Middle magnetosphere 


N e 1 cm~ 3 


0.5 


0.02 


0.01 


W th 1 keV 


0.05 


0.25 


2.5 


illio / M m ~ 2 


0.095 


0.0085 


0.013 


Efo / mWm -2 


0.0095 


0.0042 


0.067 
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Figure 1: (a) Thermospheric and plasma angular velocity profiles for the transient com- 
pression cases as a function of ionospheric latitude. Thick lines represent thermospheric 
profiles whilst thin lines represent plasma profiles. The solid lines represent case CS 
(steady state before compression) whilst the dashed and dot-dashed lines indicate cases 
CH (system at minimum disc radius) and CF (system just returned to initial disc radius) 
respectively. The magnetospheric regions (region III shaded) are labelled and separated 
by the black dotted lines. The magnetically mapped location of Io on the ionosphere is 
marked and labelled, (b) Thermosphere and plasma angular velocity profiles for the tran- 
sient expansion cases as a function of ionospheric latitude. The line styles are the same as 
(a) but the cases are now ES, EH and EF respectively, where 'E' denotes expansion, and 
the 'S', 'H' and 'F' symbols represent the same phases of the event as for Fig. [lp. This 
Figure is the same as Fig. 2 in Paper 1. 
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Figure 2: (a) FAC densities in the high latitude ionospheric region for our transient 
compression cases. The solid line represents case CS whilst the dashed and dot-dashed 
lines indicate cases CH and CF respectively. The conjugate magnetospheric regions (region 
III is shaded) are separated by dotted black lines and labelled. The latitudinal size of an 
ACS-SBC HST pixel located near the main auroral emission is represented by the dark 
grey box. This Figure is the same as Fig. 4 in Paper 1. (b) Shows the latitudinal variation 
of the precipitating electron flux for the transient compression cases. The colour codes and 
in plot labels are the same as (a) . The solid grey line indicates the limit of detectability 
of the HST flCowley et aHp007|). 
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Figure 3: (a)-(c) show the variation of thermospheric azimuthal velocity (colour scale) 
in the corotating reference frame for cases CS, CH and CF respectively (left to right). 
Positive values (dark red) indicate super-corotation, whilst negative values (light red to 
blue) indicate sub-corotation. The arrows show the direction of meridional flow and the 
white line indicates the locus of rigid corotation. The solid black encloses regions of super- 
corotation (> 25ms -1 ) and the dashed white line encloses regions that are sub-corotating 
at a rate > — 2500 m s -1 . The magnetospheric regions (region III shaded) are separated by 
dotted black lines and labelled, (d)-(f) show the meridional velocity in the thermospherc 
for cases CS-CF. The colour scale indicates the speed of flows. All other labels and symbols 
are as for (a)-(c). (g) shows the thermospheric temperature distributions for case CS whilst 
(h)-(i) show the temperature difference between cases CH and CS; and between cases CF 
and CS. Other labels are as in (a)-(c). This Figure is the same as Fig. 6 in Paper 1. 



31 




Figure 4: (a)-(c) shows the variation of atmospheric heating terms with altitude, latitude 
and temperature (colour bar) for cases CS, CH and CF (left to right). The contours 
enclose regions where heating/kinetic energy rates exceed 20Wkg~ 1 . Ion drag, Joule 
heating, vertical and horizontal advection of energy, adiabatic heating/cooling, viscous 
heating and heat conduction (vertical and turbulent) are represented by blue, red, yellow 
and magenta, green, cyan and orange lines. The magnetospheric regions are separated 
and labelled, (d)-(f) show the variation of atmospheric cooling terms where the contours 
enclose regions where heating/kinetic energy are decreasing (cooling) with rates exceeding 
20Wkg _1 . All colours and labels are as in (a)-(c). (g)-(i) show how the power per 
unit area varies for our transient compression cases. The blue line represents total power 
which is the sum of magnetospheric power (red line) and atmospheric power (green line) ; 
atmospheric power is the sum of both Joule heating (black solid line) and ion drag (cyan 
solid line). Other labels are as for (a)-(c). 
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Figure 5: Integrated ionospheric powers per hemisphere for cases CS-CF are represented 
in this figure. Ion drag is represented by blue bars, Joule heating by green bars, magneto- 
spheric power by red bars and total (sum of all above) is represented by orange bars. See 
text for further detail. 
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Figure 6: FAC densities in the high latitude region for our transient expansion cases. The 
solid line represents case ES whilst the dashed and dot-dashed lines indicate cases EH and 
EF respectively. The magnetospheric regions (region III is shaded) are separated by dotted 
black lines and labelled. The latitudinal size of an ACS-SBC HST pixel is indicated by 
the dark grey box. This Figure is the same as Fig. 8 in Paper 1. (b) Shows the latitudinal 
variation of the precipitating electron flux for the transient expansion cases. The colour 
codes and in plot labels are the same as (a). 



33 



Azimuthal Flow for ES ms" 1 Azimuthal Flow for EH ms" 1 Azimuthal Flow tot EF 




65 70 75 BO 85 90 65 70 75 80 85 90 65 70 75 80 85 90 

Latitude / ° Latitude / ° Latitude / ° 



Figure 7: (a)-(c) show the variation of thermospheric azimuthal velocity (colour scale) in 
the corotating reference frame for cases ES-EF respectively (left to right). Positive values 
(dark red) indicate super-corotation, whilst negative values (light red to blue) indicate 
sub-corotation. The arrows show the direction of meridional flow and the white line 
indicates the locus of rigid corotation. The solid black encloses regions of super-corotation 
(> 25 ms -1 ) and the dashed white line encloses regions that are sub-corotating at a rate 
> —1750 ms -1 . The magnetospheric regions (region III is shaded) are separated by dotted 
black lines and labelled, (d)-(f) show the meridional velocity in the thermosphere for cases 
ES-EF. The colour scale indicates the speed of flows. All other labels and are as for (a)-(c). 
(g) shows the thermospheric temperature distributions for case ES whilst (h)-(i) show the 
temperature difference between cases EH and EF and case ES. All labels are as in (a)-(c). 
This Figure is the same as Fig. 9 in Paper 1. 
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Figure 8: (a)-(c) shows the variation of atmospheric heating terms with altitude, latitude 
and temperature (colour bar) for cases ES, EH and EF (left to right). The contours 
enclose regions where heating/kinetic energy rates exceed 20Wkg~ 1 . Ion drag, Joule 
heating, vertical and horizontal advection of energy, adiabatic heating/cooling, viscous 
heating and heat conduction (vertical and turbulent) are represented by blue, red, yellow 
and magenta, green, cyan and orange lines. The magnetospheric regions are separated 
and labelled, (d)-(f) show the variation of atmospheric cooling terms where the contours 
enclose regions where heating/kinetic energy are decreasing (cooling) with rates exceeding 
20Wkg _1 . All colours and labels are as in (a)-(c). (g)-(i) show how the power per 
unit area varies for our transient expansion cases. The blue line represents total power 
which is the sum of magnetospheric power (red line) and atmospheric power (green line) ; 
atmospheric power is the sum of both Joule heating (black solid line) and ion drag (cyan 
solid line). All other labels are as for (a)-(c). 
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Figure 9: Integrated ionospheric powers per hemisphere for cases ES-EF are represented 
in this figure. Ion drag is represented by blue bars, Joule heating by green bars, magne- 
tospheric power by red bars and total (sum of all above) is represented by orange bars. 
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